Supersymmetric manipulation of quasienergy states. Application to the Berry phase 
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Time-dependent supersymmetry allows one to delete quasienergy levels for time-periodic Hamil- 
tonians and to create new ones. We illustrate this by examining an exactly solvable model related 
to the simple harmonic oscillator with a time-varying frequency. For an interesting nonharmonic 
example we present the change of the Berry phase due to a supersymmetry transformation. 
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One may say without exaggeration that non- 
relativistic quantum mechanics is essential to under- 
standing the basic laws of nature. Thus, the Schrodinger 
equation is one of the most fundamental constructs of 
modern natural science and over the past century, as 
evidenced by a vast literature, great progress has been 
made in developing various consequences of the time- 
independent Schrodinger equation describing stationary 
processes. Key recent advances concern the manipulation 
of bound states 0. In the real world, however, nearly 
everything changes with time, and stationary processes 
are exceptional (perhaps even impossible), but there is 
ample evidence for believing that there will soon be sim- 
ilar advances in the study of non-stationary processes 
and the consequences of the time- dependent Schrodinger 
equation. The end of the last century saw some progress 
in this area: an extensive investigation into the proper- 
tics of time periodic Hamiltonians, their dynamical in- 
variants Hand symmetries ^ dynamical pulsed coher- 
ent states y , Floquet quanta 5] and application in optics 
0. An area that is now attracting interest is calculat- 
ing the Berry, or geometric, phase for time dependent 
Hamiltonians, although so far the only system treated in 
detail is the, exactly solvable, harmonic oscillator with 
time dependent mass and/or frequency 0- 

It has long been known how to construct new exactly 
solvable time-independent Hamiltonians from an initial 
one by means of the Darboux, or supersymmetry, trans- 
formation (see e.g. H,|9]), and the procedure has recently 
been generalized to include time dependence Al- 
though stationary state energies do not exist for a time- 
dependent Hamiltonian, if the time dependence is peri- 
odic, quasienergies play a similar role [ij, 0| . In this 
letter, based on the simple harmonic oscillator model, 
having a time-periodic frequency, we show that this anal- 
ogy is much deeper than previously recognized and that 
it can be extended to the level of supersymmetry. In par- 



ticular, we show that processes such as the creation and 
annihilation of quasienergies may be effected easily by 
time-dependent supersymmetry (or Darboux) transfor- 
mations. This also opens the door for performing more 
complex operations on quasienergy states by applying 
integral transformations which are essentially inverses to 
the differential supersymmetry transformations 0, 0,0] ■ 

Since exact eigenstates are known for these new time- 
dependent Hamiltonians their remarkable properties, 
such as the the Berry phase, can be investigated in detail. 
We demonstrate here that the Berry phase, for systems 
generated from our model harmonic oscillator Hamilto- 
nian, simply acquires an additive correction under a Dar- 
boux transformation. 

We begin with the Schrodinger equation for the Hamil- 
tonian 



ho^ -dl+uj'^{t)x'^ , uj(t + T)=uj{t) 



(1) 



where dx = d/dx and appropriate units have been chosen. 
A complete set of solutions normalized to unity on the 
real Hne is [H 

iPJx. t) - exp {z^z') /„(x, t) (2) 

where 

fnix,t) = Nn-/-^/^H„{z)e-'"/\ z{x,t) = x/^/8^, 

(3) 

se, Hn{z) is a Her- 



Here 7V„ = {T^+^n\^/2^)-^/^ , 
mite polynomial and e(t), e{t) are two linearly indepen- 
dent solutions of the classical equation of motion for the 
harmonic oscillator 
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e{t) + Auj'^{t)e{t) = 0. 



(4) 



A dot over a symbol represents the derivative with re- 
spect to time and we choose these solutions so that their 
Wronskian i{t)lit) — e{t)e{t) — ^. In general, the solu- 
tions of l@J are of three types: (i) both e and e are stable, 
(ii) they are both unstable, (iii) one is stable the other is 
unstable. It is known that only in the first case the 
functions have the property 
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ipnix, t + T) = e 



■>Pnix,t) , 



(5) 
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where En = (n + is called the quasienergy and 5 
is defined by the equation 



e{t + T) = e{t)e 



IST 



(6) 



According to Floquet's theorem a (complex) solution e{t) 
to equation satisfying © can always be found in case 
(i). We choose the complex conjugate of e to be the 
second, linearly independent, solution to Here and 
below a bar over a quantity denotes its complex conjugate 
and we note that ^{t) is real. It is clear from Q that 
quasienergies are defined modulo 2t: and if they are 
comensurable with the functions ■i/'n are periodic. 

A further useful property of the functions Q is that 
there exist ladder operators 



edx — isx/2, 



-edx + iexl2 



(7) 



such that 



These functions are eigenfunctions of the symmetry oper- 
ator g = [aa^+a^a)/2^ gipn = gnipn with gn = |(2n-|-l). 

Next we shall apply the approach presented in for 
constructing new exactly solvable time-dependent Hamil- 
tonians having the same quasienergy spectrum, with the 
possible exception of a few levels, as a given one. The 
procedure is simply the time-dependent generalization 
of the usual supersymmetry construction A new 

hamiltonian hi is defined by means of a nodeless solu- 
tion u{x,t), called the transformation function, to the 
initial Schrodinger equation, idtu{x,t) ~ hQu(x,t), sub- 
ject to the additional condition d^{\iiu/u) — 0, which 
guarantees that hi has the form 



hi ^ hQ-dl\n\u{x,t)\^ 



(9) 



Solutions of the corresponding Schrodinger equation are 
obtained by applying the differential operator 



L = Li{t)[^dx+d\nu{x,t)/dx] 



to 



ipn{x,t) ^ Mn^^^LTPn{x,t) 



(10) 



(11) 



where the time-independent factor A/„ guarantees 
that the functions are normalized to unity. The op- 
erator L is defined in terms of the same function u{x,t) 
as in (9) and 



Liit) = cxp 



2 y Im {dllnu) dt 



(12) 



The operator i+ = Li{t)[dx + d\nu{x,t) / dx], adjoint to 
L, realizes the transformation in the opposite direction, 
from solutions of the Schrodinger equation with Hamil- 
tonian hi to those of the Hamiltonian /iq. Thus, the 



superposition L+L is a symmetry operator for the initial 
Schrodinger equation and the function u{x, t) is an eigen- 
function of this operator: L^Lu[x,t) = guu{x,t). The 
normalization factor in is equal to its mean value 
M„ = (i^„|L+L|V'„). 

Among the functions |(2Jl only is nodeless and suit- 
able for using as a transformation function; it produces 
only a shift of the Hamiltonian by an x-independent 
value. This is simply a manifestation of the well-known 
shape-invariance property |Q in the time-dependent case. 
Any other function Uk{x, t) = ipk{x,t), k > taken as the 
transformation function will produce a potential with k 
poles corresponding to the zeros of ipk{x, t) which clearly 
has no physical meaning if the variable x runs over the 
whole real line. The transformed Hamiltonian hi can be 
taken as the initial one for the next transformation step 
and if this is realized with the transformation function 
Uk+i{x,t) — Luk+i{x,t) all the poles are removed and 



the resulting Hamiltonian h^2^'' = ha — A2''' {x, t) is phys- 
''Pn+i ■ (8) ically admissible. For the potential difference one gets 



Af\x,t) 



1 
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Jk{z) \Jk{z) 



(13) 



Here 



and for the first three of these functions one has the sim- 
ple expressions: 

Jo(z) = l, Ji(z) = 2z2 + 1, J2{z)^Az^ + i. 

The Hamiltonian ft,2 has the same system of quasiener- 
gies as ha except for those corresponding to rt = A; and 
n = A; + 1 which are now deleted. We can repeat this 
process without any restrictions and get a Hamiltonian 
having an initial quasienergy spectrum with any number 
of lacunae composed of two adjacent levels. The poten- 



tial = ^2(^)^2 



- A^^^ looks like a harmonic oscillator 



potential, at the bottom of which there are k additional 
minima. The behavior of the k — 2 case, which is typical, 
is sketched in Fig. 1 together with the harmonic oscillator 
potential. 

The opposite process, the creation of new quasienergy 
levels, is also possible. For this purpose we need unphysi- 
cal solutions of the initial Schrodinger equation, which do 
not belong to the Hilbert space and have growing asymp- 
totic behavior as |a;| ^ 00 at any fixed moment of time. 
It is not difficult to see that the functions 

Uk{x,t)=^-'/*ie/lf/'+'^'Hk{iz)e^^^''+'/^^^' (14) 

have this property. For any even value k = 21 these 
are nodeless and suitable for using as transformation 
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FIG. 1: Potentials V^^\x,t) at t = (curve 2a) and at t ^ 
T/2 (curve 2b) together with the harmonic oscillator potential 
(curves la and lb respectively). 
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FIG. 2: Potential V-^^\x,t) at i = (curve 2a) and at t = 
T/2 (curve 2b) together with the harmonic oscillator potential 
(curves la and lb respectively). 



functions. They produce the new Hamiltonians h\ 
—9^ + Vi^''\x, t), / = 0, 1, . . . with potentials 



(20 



where 



,(20/ 



A 



(20 
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1 + 4/(2^-1)^^^- 



q2i-i{z) 



(16) 

and qk{z) = {—1)^2 ^/'^Hk{iz). Using the recursion rela- 
tion for Hermite polynomials one finds 

qniz) ^ 1, qi{z) = \/2z, qk+i{z) = \/2zqk{z) + kqk-i{z) . 

The fact that a new quasienergy level is created by 
this process follows from the property that the function 
V — l/{LiUk) is a square integrable solution of the trans- 
formed Schrodinger equation. It is easy to see that it 
corresponds to the quasienergy E = —6{k + which, 
in general, is different from all the other quasienergies 
En = S{n + i), n = 0, 1, . . .. The case fc = repro- 
duces the harmonic oscillator Hamiltonian shifted by an 
x-independent quantity. The first nontrivial case cor- 
responds to k = 2. We display a typical potential, 
Vl^\x,t), in Fig. 2. 

After applying the operator (|l()|l to the functions (O 
one gets solutions of the transformed equation corre- 
sponding to the quasienergies En- For the case fc = 2 



they can easily be expressed in terms of '0„(a;,i): 



iPn{x,t) 



Vn+3 



[V^rTT(e/e)i/2V„+i(a;,i) 



(17) 



The normalization factor here is calculated by noting that 
the symmetry operator L^L is simply the shifted g op- 
erator, L+L = 3-1-5/8. 

We now turn to the transformation of the Berry phase 
and consider for simplicity the first nontrivial case k — 2. 
In these calculations we are using the standard approach 
(see e.g. 0). The phase change of a solution of the 
Schrodinger equation 'il'{x, t) during time T is given in 
terms of the mean energy value 



X = 



(18) 



and the Berry phase /3 is obtained from this overall phase 
change, to be extracted from (2) for the harmonic oscil- 
lator potential and from (17) for transformed potentials, 
by subtracting the dynamical component. 
For the functions in (01 one gets 

i d 1 d"^ 

ii'nlti'n) = -A2n + l)-ln(e/e) - ^(V-nla^'lV-n) n^lnT 



dt 
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df^ 



(19) 

where the mean value of the square coordinate is 

(V-nk^lV-n) =47(2n+l). 

The structure of the functions tfn{x,t) is similar to 
those in (0) and therefore the average energy is given by 
()19|l with the replacement ipn fn- To calculate the 
mean value of the square coordinate for the functions 
(|17|l we use the property L+L = g -I- 5/8 and the fact 
that the functions Q are eigenfunctions for g to get 



4^/7 



The first integral here is standard and the second has the 
expression (ipnlxltpn) ~ -I- 3 — /„(^)], where 



Hl{x) 



dx 



which obeys the recursion relation 

In{a) = -2In-i{a) + 4(n - lfln-2{a) - 4a/;_i(a), 

which follows directly from that for the Hermite polyno- 
mials. The initial values are /o(a) = -^e^erfcy^ and 

/i(a) — A^/tt ~ AaIo{a). Therefore, the mean energy in 
the states (fTTjl is given by 



d^ 



After being integrated over a period T the first term in 
(|19|l gives us the overall phase change both for the initial 
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states and for the transformed ones (see l)20|l ): the Berry 
phase is determined only by the second term. Thus, the 
Berry phases /3° for all states tpn{x, t) are determined by 
the Berry phase for the ground state: /3JJ = (2n + 1)(3q, 



— dt. 



For the transformed Berry phase from H20|l one gets (3}^ — 
/3° + 2[1 - '-^0^]I3^. Hence, once the quantity /3g is 
known, we can easily calculate the Berry phase both for 

(2) 

the harmonic oscillator and for the Hamiltonian h\ . 

As a numerical illustration we have chosen a model 
for which analytic solutions of equation Q are available: 

Lj(t) — y^o^Q — ^p{t + Wi). Here ujq is a parameter of the 
model along with the real (a;,- = T) and imaginary (wi) 
periods of the Weierstrass p function. If we eliminate 
the parameter ojq in favor of d given by p{d) = — 4a;g, 
the solutions of the equation Q have the form 



e{t) 



d) 



cr(t + LU,) 



Here a and ( are (non-elliptic) Weierstrass functions. 
Figures 1 and 2 are plotted with ujr ~ —ii^i ~ 2 and 
wo = 0.5978. We have found also the value /?§ = -0.0149. 

We make one further comment concerning the time- 
dependent supersymmetry underlying our approach. The 
transformation operators are related to eigenfunctions 
(not necessarily "physical") of the symmetry operator 
g = \{a^a -f aa+), guk — guU^, gu — 5/8, which ac- 
counts for the factorization L+L = g — g^- When the or- 
der of transformation operators is interchanged, one gets 
a symmetry operator for the transformed Schrodinger 
equation g, LL^ = g — gu- The operators g and g 
are supersymmetric partners from which a matrix op- 
erator Q = diag(g,g) can be constructed. It acts in 
the space spanned by the basis vectors '^IP — {ipn,OY 
and "^n^ = (0, ifnY where the superscript "t" stands for 
transposition and the functions 4'n and ipn are given in 
(0) and (|17|l respectively. So, just as in the usual su- 
persymmetric approach, one has a two-fold degenerate 
spectrum except for the ground state level which is non- 
degenerate. This means that we have constructed here 
a model with unbroken supersymmetry. Finally we note 
that supercharge operators, closing superalgebra, can be 
constructed as usual with the help of the transformation 
operators L and L+. 

In summary, in this note we have constructed, in prin- 
ciple, an infinite number of time-periodic Hamiltonians 
having almost unlimited complexity, for which the Berry 
phase can be determined explicitly and have illustrated 
the procedure for a Hamiltonian whose time dependence 
is given by the square root of an elliptic function. This 



vastly extends the set of previously known cases, all based 
on the simple Harmonic oscillator. We feel that this 
opens the way for a systematic investigation of the Berry 
phase for one dimensional quantum systems. 
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